Exome-wide variation in a diverse barley panel reveals genetic associations with ten agronomic traits in Eastern landraces.

Barley (Hordeum vulgare ssp. vulgare) is one of the first crops to be domesticated and is adapted to a wide range of environments. Worldwide barley germplasm collections possess valuable allelic variations that could further improve barley productivity. Although barley genomics has offered a global picture of allelic variation among varieties and its association with various agronomic traits, polymorphisms from East Asian varieties remain scarce. In this study, we analyze exome polymorphisms in a panel of 274 barley varieties collected worldwide, including 137 varieties from East Asian countries and Ethiopia. We reveal the underlying population structure and conduct genome-wide association studies for 10 agronomic traits. Moreover, we examin genome-wide associations for traits related to grain size such as awn length and glume length. Our results demonstrate the value of diverse barley germplasm panels containing Eastern varieties, highlighting their distinct genomic signatures relative to Western subpopulations.

Transcriptome Analysis of Chloris virgata, Which Shows the Fastest Germination and Growth in the Major Mongolian Grassland Plant.

Plants in Mongolian grasslands are exposed to short, dry summers and long, cold winters. These plants should be prepared for fast germination and growth activity in response to the limited summer rainfall. The wild plant species adapted to the Mongolian grassland environment may allow us to explore useful genes, as a source of unique genetic codes for crop improvement. Here, we identified the Chloris virgata Dornogovi accession as the fastest germinating plant in major Mongolian grassland plants. It germinated just 5 h after treatment for germination initiation and showed rapid growth, especially in its early and young development stages. This indicates its high growth potential compared to grass crops such as rice and wheat. By assessing growth recovery after animal bite treatment (mimicked by cutting the leaves with scissors), we found that C. virgata could rapidly regenerate leaves after being damaged, suggesting high regeneration potential against grazing. To analyze the regulatory mechanism involved in the high growth potential of C. virgata, we performed RNA-seq-based transcriptome analysis and illustrated a comprehensive gene expression map of the species. Through de novo transcriptome assembly with the RNA-seq reads from whole organ samples of C. virgata at the germination stage (2 days after germination, DAG), early young development stage (8 DAG), young development stage (17 DAG), and adult development stage (28 DAG), we identified 21,589 unified transcripts (contigs) and found that 19,346 and 18,156 protein-coding transcripts were homologous to those in rice and Arabidopsis, respectively. The best-aligned sequences were annotated with gene ontology groups. When comparing the transcriptomes across developmental stages, we found an over-representation of genes involved in growth regulation in the early development stage in C. virgata. Plant development is tightly regulated by phytohormones such as brassinosteroids, gibberellic acid, abscisic acid, and strigolactones. Moreover, our transcriptome map demonstrated the expression profiles of orthologs involved in the biosynthesis of these phytohormones and their signaling networks. We discuss the possibility that C. virgata phytohormone signaling and biosynthesis genes regulate early germination and growth advantages. Comprehensive transcriptome information will provide a useful resource for gene discovery and facilitate a deeper understanding of the diversity of the regulatory systems that have evolved in C. virgata while adapting to severe environmental conditions.

Cellular and transcriptomic analyses reveal two-staged chloroplast biogenesis underpinning photosynthesis build-up in the wheat leaf.

BACKGROUND: The developmental gradient in monocot leaves has been exploited to uncover leaf developmental gene expression programs and chloroplast biogenesis processes. However, the relationship between the two is barely understood, which limits the value of transcriptome data to understand the process of chloroplast development. RESULTS: Taking advantage of the developmental gradient in the bread wheat leaf, we provide a simultaneous quantitative analysis for the development of mesophyll cells and of chloroplasts as a cellular compartment. This allows us to generate the first biologically-informed gene expression map of this leaf, with the entire developmental gradient from meristematic to fully differentiated cells captured. We show that the first phase of plastid development begins with organelle proliferation, which extends well beyond cell proliferation, and continues with the establishment and then the build-up of the plastid genetic machinery. The second phase is marked by the development of photosynthetic chloroplasts which occupy the available cellular space. Using a network reconstruction algorithm, we predict that known chloroplast gene expression regulators are differentially involved across those developmental stages. CONCLUSIONS: Our analysis generates both the first wheat leaf transcriptional map and one of the most comprehensive descriptions to date of the developmental history of chloroplasts in higher plants. It reveals functionally distinct plastid and chloroplast development stages, identifies processes occurring in each of them, and highlights our very limited knowledge of the earliest drivers of plastid biogenesis, while providing a basis for their future identification.

Phosphate or nitrate imbalance induces stronger molecular responses than combined nutrient deprivation in roots and leaves of chickpea plants.

The negative effects of phosphate (Pi) and/or nitrate (NO3- ) fertilizers on the environment have raised an urgent need to develop crop varieties with higher Pi and/or nitrogen use efficiencies for cultivation in low-fertility soils. Achieving this goal depends upon research that focuses on the identification of genes involved in plant responses to Pi and/or NO3- starvation. Although plant responses to individual deficiency in either Pi (-Pi/+NO3- ) or NO3- (+Pi/-NO3- ) have been separately studied, our understanding of plant responses to combined Pi and NO3- deficiency (-Pi/-NO3- ) is still very limited. Using RNA-sequencing approach, transcriptome changes in the roots and leaves of chickpea cultivated under -Pi/+NO3- , +Pi/-NO3- or -Pi/-NO3- conditions were investigated in a comparative manner. -Pi/-NO3- treatment displayed lesser effect on expression changes of genes related to Pi or NO3- transport, signalling networks, lipid remodelling, nitrogen and Pi scavenging/remobilization/recycling, carbon metabolism and hormone metabolism than -Pi/+NO3- or +Pi/-NO3- treatments. Therefore, the plant response to -Pi/-NO3- is not simply an additive result of plant responses to -Pi/+NO3- and +Pi/-NO3- treatments. Our results indicate that nutrient imbalance is a stronger stimulus for molecular reprogramming than an overall deficiency.

Identification of effector candidate genes of Rhizoctonia solani AG-1 IA expressed during infection in Brachypodium distachyon.

Rhizoctonia solani is a necrotrophic phytopathogen belonging to basidiomycetes. It causes rice sheath blight which inflicts serious damage in rice production. The infection strategy of this pathogen remains unclear. We previously demonstrated that salicylic acid-induced immunity could block R. solani AG-1 IA infection in both rice and Brachypodium distachyon. R. solani may undergo biotrophic process using effector proteins to suppress host immunity before necrotrophic stage. To identify pathogen genes expressed at the early infection process, here we developed an inoculation method using B. distachyon which enables to sample an increased amount of semi-synchronous infection hyphae. Sixty-one R. solani secretory effector-like protein genes (RsSEPGs) were identified using in silico approach with the publicly available gene annotation of R. solani AG-1 IA genome and our RNA-sequencing results obtained from hyphae grown on agar medium. Expression of RsSEPGs was analyzed at 6, 10, 16, 24, and 32 h after inoculation by a quantitative reverse transcription-polymerase chain reaction and 52 genes could be detected at least on a single time point tested. Their expressions showed phase-specific patterns which were classified into 6 clusters. The 23 RsSEPGs in the cluster 1-3 and 29 RsSEPGs in the cluster 4-6 are expected to be involved in biotrophic and necrotrophic interactions, respectively.

BdWRKY38 is required for the incompatible interaction of Brachypodium distachyon with the necrotrophic fungus Rhizoctonia solani.

Rhizoctonia solani is a soil-borne necrotrophic fungus that causes sheath blight in grasses. The basal resistance of compatible interactions between R. solani and rice is known to be modulated by some WRKY transcription factors (TFs). However, genes and defense responses involved in incompatible interaction with R. solani remain unexplored, because no such interactions are known in any host plants. Recently, we demonstrated that Bd3-1, an accession of the model grass Brachypodium distachyon, is resistant to R. solani and, upon inoculation with the fungus, undergoes rapid induction of genes responsive to the phytohormone salicylic acid (SA) that encode the WRKY TFs BdWRKY38 and BdWRKY44. Here, we show that endogenous SA and these WRKY TFs positively regulate this accession-specific R. solani resistance. In contrast to a susceptible accession (Bd21), the infection process in the resistant accessions Bd3-1 and Tek-3 was suppressed at early stages before the development of fungal biomass and infection machinery. A comparative transcriptome analysis during pathogen infection revealed that putative WRKY-dependent defense genes were induced faster in the resistant accessions than in Bd21. A gene regulatory network (GRN) analysis based on the transcriptome dataset demonstrated that BdWRKY38 was a GRN hub connected to many target genes specifically in resistant accessions, whereas BdWRKY44 was shared in the GRNs of all three accessions. Moreover, overexpression of BdWRKY38 increased R. solani resistance in Bd21. Our findings demonstrate that these resistant accessions can activate an incompatible host response to R. solani, and BdWRKY38 regulates this response by mediating SA signaling.

De novo transcriptome analysis reveals an unperturbed transcriptome under high cadmium conditions in the Cd-hypertolerant fern Athyrium yokoscense.

Athyrium yokoscense shows strong tolerance to cadmium exposure, even at levels that are many times greater than the toxic levels in ordinary plants. To determine the mechanism of Cd tolerance in A. yokoscense, we grew these plants under high Cd conditions and observed the tissue-specific accumulation of Cd and generation of reactive oxygen species, which is one of the major physiological responses to Cd stress. Fuchsin staining indicated the existence of a casparian strip in A. yokoscense roots, which may participate in Cd hypertolerance in A. yokoscense. Moreover, we performed RNA-seq of RNA samples from A. yokoscense plants treated with or without Cd exposure and obtained comprehensive RNA sequences as well as the Cd-responsive expression patterns of individual genes. Through de novo transcriptome assembly and gene expression analysis, we found that A. yokoscense showed normal features with no significant change in the expression levels of any transporter genes, even under high Cd exposure conditions. Our results demonstrate that A. yokoscense has an unusual mechanism that allows the invading Cd to partition into the distal roots, thus avoiding translocation of Cd into the xylem.

Genetic Factors Associated with Heading Responses Revealed by Field Evaluation of 274 Barley Accessions for 20 Seasons.

Heading time is a key trait in cereals affecting the maturation period for optimal grain filling before harvest. Here, we aimed to understand the factors controlling heading time in barley (Hordeum vulgare). We characterized a set of 274 barley accessions collected worldwide by planting them for 20 seasons under different environmental conditions at the same location in Kurashiki, Japan. We examined interactions among accessions, known genetic factors, and an environmental factor to determine the factors controlling heading response. Locally adapted accessions have been selected for genetic factors that stabilize heading responses appropriate for barley cultivation, and these accessions show stable heading responses even under varying environmental conditions. We identified vernalization requirement and PPD-H1 haplotype as major stabilizing mechanisms of the heading response for regional adaptation in Kurashiki.

Parental legacy and regulatory novelty in Brachypodium diurnal transcriptomes accompanying their polyploidy.

Polyploidy is a widespread phenomenon in eukaryotes that can lead to phenotypic novelty and has important implications for evolution and diversification. The modification of phenotypes in polyploids relative to their diploid progenitors may be associated with altered gene expression. However, it is largely unknown how interactions between duplicated genes affect their diurnal expression in allopolyploid species. In this study, we explored parental legacy and hybrid novelty in the transcriptomes of an allopolyploid species and its diploid progenitors. We compared the diurnal transcriptomes of representative Brachypodium cytotypes, including the allotetraploid Brachypodium hybridum and its diploid progenitors Brachypodium distachyon and Brachypodium stacei. We also artificially induced an autotetraploid B. distachyon. We identified patterns of homoeolog expression bias (HEB) across Brachypodium cytotypes and time-dependent gain and loss of HEB in B. hybridum. Furthermore, we established that many genes with diurnal expression experienced HEB, while their expression patterns and peak times were correlated between homoeologs in B. hybridum relative to B. distachyon and B. stacei, suggesting diurnal synchronization of homoeolog expression in B. hybridum. Our findings provide insight into the parental legacy and hybrid novelty associated with polyploidy in Brachypodium, and highlight the evolutionary consequences of diurnal transcriptional regulation that accompanied allopolyploidy.

Plastic transport systems of rice for mineral elements in response to diverse soil environmental changes.

Climate change will increase frequency of drought and flooding, which threaten global crop productivity and food security. Rice (Oryza sativa) is unique in that it is able to grow in both flooded and upland conditions, which have large differences in the concentrations and chemical forms of mineral elements available to plants. To comprehensively understand the mechanisms of rice for coping with different water status, we performed ionomics and transcriptomics analysis of the roots, nodes and leaves of rice grown in flooded and upland conditions. Focusing the analysis on genes encoding proteins involved in transport functions for mineral elements, it was found that, although rice plants maintained similar levels of each element in the shoots for optimal growth, different transporters for mineral elements were utilised for nitrogen, iron, copper and zinc to deal with different soil water conditions. For example, under flooded conditions, rice roots take up nitrogen using transporters for both ammonium (OsAMT1/2) and nitrate (OsNPF2.4, OsNRT1.1A and OsNRT2.3), whereas under upland conditions, nitrogen uptake is mediated by different nitrate transporters (OsNRT1.1B and OsNRT1.5A). This study shows that rice possesses plastic transport systems for mineral elements in response to different water conditions (upland and flooding).

Genetic Variation for Seed Metabolite Levels in Brachypodium distachyon.

Metabolite composition and concentrations in seed grains are important traits of cereals. To identify the variation in the seed metabolotypes of a model grass, namely Brachypodium distachyon, we applied a widely targeted metabolome analysis to forty inbred lines of B. distachyon and examined the accumulation patterns of 183 compounds in the seeds. By comparing the metabolotypes with the population structure of these lines, we found signature metabolites that represent different accumulation patterns for each of the three B. distachyon subpopulations. Moreover, we found that thirty-seven metabolites exhibited significant differences in their accumulation between the lines Bd21 and Bd3-1. Using a recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21, we identified the quantitative trait loci (QTLs) linked with this variation in the accumulation of thirteen metabolites. Our metabolite QTL analysis illustrated that different genetic factors may presumably regulate the accumulation of 4-pyridoxate and pyridoxamine in vitamin B6 metabolism. Moreover, we found two QTLs on chromosomes 1 and 4 that affect the accumulation of an anthocyanin, chrysanthemin. These QTLs genetically interacted to regulate the accumulation of this compound. This study demonstrates the potential for metabolite QTL mapping in B. distachyon and provides new insights into the genetic dissection of metabolomic traits in temperate grasses.

Characterization of sulfur-compound metabolism underlying wax-ester fermentation in Euglena gracilis.

Euglena gracilis is a microalga, which has been used as a model organism for decades. Recent technological advances have enabled mass cultivation of this species for industrial applications such as feedstock in nutritional foods and cosmetics. E. gracilis degrades its storage polysaccharide (paramylon) under hypoxic conditions for energy acquisition by an oxygen-independent process and accumulates high amount of wax-ester as a by-product. Using this sequence of reactions referred to as wax-ester fermentation, E. gracilis is studied for its application in biofuel production. Although the wax-ester production pathway is well characterized, little is known regarding the biochemical reactions underlying the main metabolic route, especially, the existence of an unknown sulfur-compound metabolism implied by the nasty odor generation accompanying the wax-ester fermentation. In this study, we show sulfur-metabolomics of E. gracilis in aerobic and hypoxic conditions, to reveal the biochemical reactions that occur during wax-ester synthesis. Our results helped us in identifying hydrogen sulfide (H2S) as the nasty odor-producing component in wax-ester fermentation. In addition, the results indicate that glutathione and protein degrades during hypoxia, whereas cysteine, methionine, and their metabolites increase in the cells. This indicates that this shift of abundance in sulfur compounds is the cause of H2S synthesis.

Computer vision-based phenotyping for improvement of plant productivity: a machine learning perspective.

Employing computer vision to extract useful information from images and videos is becoming a key technique for identifying phenotypic changes in plants. Here, we review the emerging aspects of computer vision for automated plant phenotyping. Recent advances in image analysis empowered by machine learning-based techniques, including convolutional neural network-based modeling, have expanded their application to assist high-throughput plant phenotyping. Combinatorial use of multiple sensors to acquire various spectra has allowed us to noninvasively obtain a series of datasets, including those related to the development and physiological responses of plants throughout their life. Automated phenotyping platforms accelerate the elucidation of gene functions associated with traits in model plants under controlled conditions. Remote sensing techniques with image collection platforms, such as unmanned vehicles and tractors, are also emerging for large-scale field phenotyping for crop breeding and precision agriculture. Computer vision-based phenotyping will play significant roles in both the nowcasting and forecasting of plant traits through modeling of genotype/phenotype relationships.

Statistical and Machine Learning Approaches to Predict Gene Regulatory Networks From Transcriptome Datasets.

Statistical and machine learning (ML)-based methods have recently advanced in construction of gene regulatory network (GRNs) based on high-throughput biological datasets. GRNs underlie almost all cellular phenomena; hence, comprehensive GRN maps are essential tools to elucidate gene function, thereby facilitating the identification and prioritization of candidate genes for functional analysis. High-throughput gene expression datasets have yielded various statistical and ML-based algorithms to infer causal relationship between genes and decipher GRNs. This review summarizes the recent advancements in the computational inference of GRNs, based on large-scale transcriptome sequencing datasets of model plants and crops. We highlight strategies to select contextual genes for GRN inference, and statistical and ML-based methods for inferring GRNs based on transcriptome datasets from plants. Furthermore, we discuss the challenges and opportunities for the elucidation of GRNs based on large-scale datasets obtained from emerging transcriptomic applications, such as from population-scale, single-cell level, and life-course transcriptome analyses.

Benzothiadiazole, a plant defense inducer, negatively regulates sheath blight resistance in Brachypodium distachyon.

Plant defense inducers that mimic functions of the plant immune hormone salicylic acid (SA) often affect plant growth. Although benzothiadiazole (BTH), a synthetic analog of SA, has been widely used to protect crops from diseases by inducing plant defense responses, we recently demonstrated that SA, but not BTH, confers resistance against Rhizoctonia solani, the causal agent of sheath blight disease, in Brachypodium distachyon. Here, we demonstrated that BTH compromised the resistance of Bd3-1 and Gaz4, the two sheath blight-resistant accessions of B. distachyon, which activate SA-dependent signaling following challenge by R. solani. Moreover, upon analyzing our published RNA-seq data from B. distachyon treated with SA or BTH, we found that BTH specifically induces expression of genes related to chloroplast function and jasmonic acid (JA) signaling, suggesting that BTH attenuates R. solani resistance by perturbing growth-defense trade-offs and/or by inducing a JA response that may increase susceptibility to R. solani. Our findings demonstrated that BTH does not work as a simple mimic of SA in B. distachyon, and consequently may presumably cause unfavorable side effects through the transcriptional alteration, particularly with respect to R. solani resistance.

Gene Co-expression Network Analysis Suggests the Existence of Transcriptional Modules Containing a High Proportion of Transcriptionally Differentiated Homoeologs in Hexaploid Wheat.

Genome duplications aid in the formation of novel molecular networks through regulatory differentiation of the duplicated genes and facilitate adaptation to environmental change. Hexaploid wheat, Triticum aestivum, contains three homoeologous chromosome sets, the A-, B-, and D-subgenomes, which evolved through interspecific hybridization and subsequent whole-genome duplication. The divergent expression patterns of the homoeologs in hexaploid wheat suggest that they have undergone transcriptional and/or functional differentiation during wheat evolution. However, the distribution of transcriptionally differentiated homoeologs in gene regulatory networks and their related biological functions in hexaploid wheat are still largely unexplored. Therefore, we retrieved 727 publicly available wheat RNA-sequencing (RNA-seq) datasets from various tissues, developmental stages, and conditions, and identified 10,415 expressed homoeologous triplets. Examining the co-expression modules in the wheat transcriptome, we found that 66% of the expressed homoeologous triplets possess all three homoeologs grouped in the same co-expression modules. Among these, 15 triplets contain co-expressed homoeologs with differential expression levels between homoeoalleles across ≥ 95% of the 727 RNA-seq datasets, suggesting a consistent trend of homoeolog expression bias. In addition, we identified 2,831 differentiated homoeologs that showed gene expression patterns that deviated from those of the other two homoeologs. We found that seven co-expression modules contained a high proportion of such differentiated homoeologs, which accounted for ≥ 20% of the genes in each module. We also found that five of the co-expression modules are abundantly composed of genes involved in biological processes such as chloroplast biogenesis, RNA metabolism, putative defense response, putative posttranscriptional modification, and lipid metabolism, thereby suggesting that, the differentiated homoeologs might highly contribute to these biological functions in the gene network of hexaploid wheat.

Transcriptional alterations during proliferation and lignification in Phyllostachys nigra cells.

Highly-lignified culms of bamboo show distinctive anatomical and mechanical properties compared with the culms of other grass species. A cell culture system for Phyllostachys nigra has enabled investigating the alterations in cellular states associated with secondary cell wall formation during its proliferation and lignification in woody bamboos. To reveal transcriptional changes related to lignification in bamboo, we analyzed transcriptome in P. nigra cells treated with the synthetic auxin 2,4-dichlorophenoxyacetic acid (2,4-D) and the synthetic cytokinin benzylaminopurine (BA) by RNA-seq analysis. We found that some genes putatively involved in cell wall biogenesis and cell division were up-regulated in response to the 2,4-D treatment, and the induction of lignification by the BA treatment was correlated with up-regulation of genes involved in the shikimate pathway. We also found that genes encoding MYB transcription factors (TFs) show correlated expression patterns with those encoding cinnamyl alcohol dehydrogenase (CAD), suggesting that MYB TFs presumably regulate secondary cell wall formation in the bamboo cells. These findings suggest that cytokinin signaling may regulate lignification in P. nigra cells through coordinated transcriptional regulation and metabolic alterations. Our results have also produced a useful resource for better understanding of secondary cell wall formation in bamboo plants.

Multiplex PCR Targeted Amplicon Sequencing (MTA-Seq): Simple, Flexible, and Versatile SNP Genotyping by Highly Multiplexed PCR Amplicon Sequencing.

Next-generation sequencing (NGS) technologies have enabled genome re-sequencing for exploring genome-wide polymorphisms among individuals, as well as targeted re-sequencing for the rapid and simultaneous detection of polymorphisms in genes associated with various biological functions. Therefore, a simple and robust method for targeted re-sequencing should facilitate genotyping in a wide range of biological fields. In this study, we developed a simple, custom, targeted re-sequencing method, designated "multiplex PCR targeted amplicon sequencing (MTA-seq)," and applied it to the genotyping of the model grass Brachypodium distachyon. To assess the practical usability of MTA-seq, we applied it to the genotyping of genome-wide single-nucleotide polymorphisms (SNPs) identified in natural accessions (Bd1-1, Bd3-1, Bd21-3, Bd30-1, Koz-1, Koz-3, and Koz-4) by comparing the re-sequencing data with that of reference accession Bd21. Examination of SNP-genotyping accuracy in 443 amplicons from eight parental accessions and an F1 progeny derived by crossing of Bd21 and Bd3-1 revealed that ~95% of the SNPs were correctly called. The assessment suggested that the method provided an efficient framework for accurate and robust SNP genotyping. The method described here enables easy design of custom target SNP-marker panels in various organisms, facilitating a wide range of high-throughput genetic applications, such as genetic mapping, population analysis, molecular breeding, and genomic diagnostics.

Homoeolog-specific activation of genes for heat acclimation in the allopolyploid grass Brachypodium hybridum.

Background: Allopolyploid plants often show wider environmental tolerances than their ancestors; this is expected to be due to the merger of multiple distinct genomes with a fixed heterozygosity. The complex homoeologous gene expression could have been evolutionarily advantageous for the adaptation of allopolyploid plants. Despite multiple previous studies reporting homoeolog-specific gene expression in allopolyploid species, there are no clear examples of homoeolog-specific function in acclimation to a long-term stress condition. Results: We found that the allopolyploid grass Brachypodium hybridum and its ancestor Brachypodium stacei show long-term heat stress tolerance, unlike its other ancestor, Brachypodium distachyon. To understand the physiological traits of B. hybridum, we compared the transcriptome of the 3 Brachypodium species grown under normal and heat stress conditions. We found that the expression patterns of approximately 26% and approximately 38% of the homoeolog groups in B. hybridum changed toward nonadditive expression and nonancestral expression, respectively, under normal condition. Moreover, we found that B. distachyon showed similar expression patterns between normal and heat stress conditions, whereas B. hybridum and B. stacei significantly altered their transcriptome in response to heat after 3 days of stress exposure, and homoeologs that were inherited from B. stacei may have contributed to the transcriptional stress response to heat in B. hybridum. After 15 days of heat exposure, B. hybridum and B. stacei maintained transcriptional states similar to those under normal conditions. These results suggest that an earlier response to heat that was specific to homoeologs originating from B. stacei contributed to cellular homeostasis under long-term heat stress in B. hybridum. Conclusions: Our results provide insights into different regulatory events of the homoeo-transcriptome that are associated with stress acclimation in allopolyploid plants.